Step-edge instability during epitaxial growth of graphene from SiC(OOOl) 
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The unique electronic properties of graphene 
offer the possibility that it could replace 
silicon when microelectronics evolves to 
nanoelectronics.[l] Graphene grown epitaxi- 
ally on silicon carbide [2] is particularly attractive 
in this regard because SiC is itself a useful 
semiconductor and, by suitable manipulation of 
the growth conditions, [sj, 0? l5l] epitaxial films 
can be produced that exhibit all the transport 
properties of ideal, two-dimensional graphene 
desired for device applications. [6, 7] Never- 
theless, there is little or no understanding of 
the actual kinetics of growth, which is likely 
to be required for future process control. As 
a step in this direction, we propose a local 
heat release mechanism to explain finger-like 
structures observed when graphene is grown by 
step flow decomposition of SiC (0001). Using 
a continuum equation of motion for the shape 
evolution of a moving step, a linear stability 
analysis predicts whether a shape perturbation 
of a straight moving step grows or decays as a 
function of growth temperature, the background 
pressure of Si maintained during growth, and 
the effectiveness of an inert buffer gas to retard 
the escape of Si atoms from the crystal surface. 
The theory gives semi-quantitative agreement 
with experiment for the characteristic separation 
between fingers observed when graphene is 
grown in a low-pressure induction furnace or 
under ultrahigh vacuum conditions. 

Figure [U is an atomic force microscope image of 2-3 ML 
of graphene grown on a commercial wafer of SiC (0001) 
at ~ 1600° C in a low-pressure induction furnace. The 
original wafer was H-etched at high temperature to pro- 
duce a regular array of steps. When heated above 
the graphitization temperature, Tg, the steps bunch into 
macrosteps (which run vertically in Figure [T]) separated 
by large nearly flat terraces. A continuous and conformal 
film of graphene covers the entire sample. A step with a 
highly ramified edge appears between each pair of adja- 
cent macrosteps. A similarly ramified step morphology 
has been reported by Hupalo and co-workers for graphene 
grown on a vicinal surface of SiC (0001) under ultrahigh 
vacuum (UHV) conditions. [ll] The thermal decomposi- 
tion process that produces graphene occurs preferentially 
at steps, where Si and C atoms are least well bonded. 
Figure [SJa) shows a step on SiC (0001) covered confor- 




FIG. 1: AFM image from Ref.lslof a SiC(OOOl) surface graphi- 
tized at T ^ 1600° C in an induction furnace with a back- 
ground Si pressure P ^ 1 x 10~^ torr. The field of view is 10 
/im. The steps descend from left to right. 



mally by one complete layer of graphene and a second 
incomplete layer. [12] The diagram reflects the fact that, 
above - 1100° C, the bulk of the SiC crystal is "buffered" 
from anything that lies above (vacuum or graphene) by 
a carbon-rich 6a/3 x 6a/3 reconstruction. The structure 
of the buffer layer is not completely known, but it ap- 
parently resembles graphene itself, albeit distorted so it 
bonds strongly to the silicon layer immediately below. 
When the temperature exceeds 1200° C, the SiC beneath 
the buffer layer of the upper terrace begins to decom- 
pose. The liberated Si atoms desorb, the upper buffer 
layer transforms to graphene, and the liberated C atoms 
recrystallize at the base of the step to extend the lower 
buffer layer. Figure [2](b) shows the advance of the second 
graphene layer. 

It is important to appreciate that surface diffusion on 
the terraces plays very little role in the growth scenario 
sketched in Figure O Therefore, despite their appear- 
ance, the fingers in Figure [T] cannot arise from a Mullins- 
Sekerka type morphological instability driven by surface 
diffusion and an asymmetry in attachment/detachment 
kinetics at the steps, [l^, ll5| We suggest an alternative 
mechanism based on the fact that the crystallization of 
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FIG. 2: Graphene growth at a step on SiC (001): (a) initial 
state; (b) final state. 



free carbon atoms into new buffer layer material cre- 
ates strong cr-bonds. The accompanying release of heat 
should transiently increase the local temperature and 
induce further decomposition. This is a positive feed- 
back mechanism that promotes further decomposition at 
points along the step where it has already begun. Since 
graphene production accompanies decomposition, Fig- 
ure [3] shows how the step edge between graphene and 
the upper buffer layer may be expected to evolve as time 
goes on. Comparison with Figure [1] identifies the light- 
colored fingers as buffer layer and the darker channels be- 
tween them as graphene. To be more quantitative, we let 
h{x, t) denote the position of the SiC step edge and take 
account of the local heating effect (phenomenologically) 
by supposing that the decomposition rate depends on the 
local curvature of the front; a higher rate for concave re- 
gions and lower rate for convex regions. Therefore, if V 
is the average velocity of the step due to decomposition, 
a suitable evolution equation for h{x,t) is 



dh 
'dt 



l^ = -V-aV^+ar^-aD^. (1) 



dx^' 



The last two terms in ([T]) are familiar from extensive stud- 
ies of the effect of capillary smoothing on the morphol- 
ogy of step edges on vicinal surfaces. lH, [l^ The second 
derivative term models evaporation-condensation events 
where atoms detach from a step edge, migrate rapidly on 
the adjacent terrace, and re- attach to the step elsewhere. 
The fourth derivative term models edge diffusion events 
where atoms migrate along the step edge itself. In ([T]), 
a is the SiC lattice constant, <j = a^^/kT (7 is the SiC 
step stiffness), F = u ex.p{—Ei/kT) is the mean rate at 
which atomic species detach from a straight SiC step {v 
is an attempt frequency) and D = a'^v ex.-p{—E2/kT) is 
the edge diffusion constant. 
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FIG. 3: Onset of instability: (a) thermal decomposition of 
SiC starts at random points along the step-edges; (b) the for- 
mation of new buffer layer releases heat, which locally acceler- 
ates further decomposition and graphitization. The topmost 
complete graphene layer is not shown. 



Kinetic theory relates the average step velocity to the 
difference between the flux of Si atoms subliming from 
the surface (as measured by the equilibrium vapor pres- 
sure of Si over SiC) and the background pressure P in 
the growth chamber: 



v = p 



(^q - P)- 



(2) 



In this expression, m is the mass of a Si atom, 
logioPeq(Pa) = 12.74 - 2.66 x 10^/T{° K),[l3 and p is 
the evaporation coefficient. [18] We note that increasing 
P shifts Tg to higher temperature and no decomposition 
occurs if P > Peq-[l3] It is straightforward to perform a 
linear stability analysis of ([1]) assuming a solution of the 
form h{x^t) = —Vt + e(t) sin(27rx/A). The perturbation 
amplitude e{t) growth (decays) exponentially when the 
wavelength A is greater (less) than a critical wavelength. 
The most unstable (fasting growing) wavelength is 
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(3) 



The constants that enter Am that require special discus- 
sion are 7, Ei, £^2, and since none of these is truly 
known for SiC. For the step stiffness of a three-bilayer 
step, we chose 7 = 1 eV/A, which is approximately 
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FIG. 4: Dependence of the fastest growing wavelength A m on 
the growth temperature T for several values of the background 
Si pressure P. 



three times the stiffness of steps on Si(lll).[16| The de- 
tachment and surface diffusion barriers were chosen as 
£^1 = 5.6 eV and E2 = 5.0 eV based on values for surface 
diffusion barriers for C atoms and Si atoms on several SiC 
surfaces given by Fissel.[20] The evaporation coefficient 
was fixed at = 0.046 to make the calculated growth 
rate equal the experimental growth rate associated with 
Figure [H 

Figure [His a plot of Am as a function of growth temper- 
ature for several values of Si background pressure. There 
is no graphene growth when T < Tq and step edges are 
absolutely stable (no fingers) when T > Ts- Between Tq 
and T^, there is a temperature window where the step 
edge is unstable. This window is largest for UHV growth 
and decreases as P increases. There is a critical pres- 
sure, P^, above which graphene growth occurs with no 
morphological instability of the SiC step edges. Below 
Ps and just above T^, there is a narrow region of sta- 
ble growth. All these features arise from the competition 
between the destabilizing effect of local heating and the 
stabilizing effect of step edge evaporation/condensation 
and edge diffusion. The capillary forces of step smoothing 
always prevail when T > Ts- The narrow temperature 
range of edge stability just above Tq occurs because the 
destabilizing curvature term in ([T]) is proportional to V, 
which is very small near Tq- Otherwise, local heating 
dominates and the scenario sketched in Figure [3] occurs. 

For the growth conditions stated in Figure [U Fig- 
ure [41(c) predicts a most unstable wavelength of ~ 0.1 
jiia. This agrees well with the separation between fin- 
gers in the AFM image. Additional experiments where 
graphene was grown in a furnace with a higher value 
of background pressure produced only straight steps, [ij 
in qualitative agreement with Figure [H^d) . Figure [4|^ a) 
corresponds to an ultrahigh vacuum environment and so 
should be relevant to the graphene growth experiments 



reported in Ref. [UJ. M T ^ 1200° C, step edge fin- 
gers were observed with a mean separation of A^, ~ 0.01 
/im, which accords well with our theory. The fingers dis- 
appeared when the UHV experiments were repeated at 
T 1350° C, [21,] which agrees only semi-quantitatively 
with Figure [H^a). On the other hand, the value of the 
evaporation coefficient /3 in ([2j) is sample-dependent and 
probably differs between the induction furnace experi- 
ment and the UHV experiment. [22] A reduction of /3 by 
only a factor of two is sufficient to close the window in 
Figure [H^a) enough to match the UHV data. Another 
process variable used in recent graphene growth exper- 
iments is the pressure of an inert _gas intentionally in- 
troduced into the growth chamber. 0, IH The presence of 
such a gas was found to tremendously improve the qual- 
ity of the graphene produced epitaxially on SiC (0001). 
Presumably, the effect of the inert gas is to retard the 
escape of Si atoms from the decomposing crystal. The 
same thing happens in the present model if, as above, we 
reduce the value of the evaporation coefficient (3 in 
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FIG. 5: Dependence of the fastest growing wavelength Am on 
the evaporation coefficient /3 for P = and T = 1400 °C. 



Accordingly, Figure [5] shows the dependence of the most 
unstable wavelength on at T = 1400° C and P = 0. 
There is a critical value, Ps^ below which the theory pre- 
dicts graphene growth without a step edge instability. 
We encourage experiments to test this prediction. 

In summary, we have shown that finger-like step edges 
observed during the growth of epitaxial graphene by 
step flow decomposition of SiC(OOOl) can be understood 
from a competition between capillary smoothing and a 
curvature-driven mechanism for step edge roughening. 
We suggest (but cannot prove) that the origin of the 
roughening is accelerated step decomposition by local 
heating at points where it has already begun. The heat 
source is the formation energy of a graphene-like buffer 
layer that accompanies the growth of graphene on this 
surface. A linear stability analysis of a step equation of 
motion that takes account of all these effects predicts the 
separation between fingers as a function of the growth 
process variables: temperature, background Si pressure, 
and the influence of an inert gas in the growth cham- 
ber. The theory agrees semi-quantitatively with growth 
experiments conducted in both a low-pressure induction 
furnace and under ultrahigh vacuum conditions. 
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